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/\ I The values and energy dependence for the survival probability <| S p> of large rapidity 

d • gaps (LRG) are calculated in a three channel model. This model includes single and double 

diffractive production, as well as elastic rescattering. It is shown that <| 5* p> decreases 
with increasing energy, in line with recent results for LRG dijet production at the Tevatron. 
This is in spite of the weak dependence on energy of the ratio {aei + o'sD)/o'tot- 
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1 Introduction 

A large rapidity gap ( LRG ) process is defined as one where no hadrons are produced in 
a sufficiently large rapidity region. Historically, both Dokshitzer et al.[0] and Bjorken[0], 
suggested utilizing LRG as a signature for Higgs production in a W-W fusion process, in 
hadron-hadron collisions. It turns out that the LRG processes give a unique opportunity 
to measure the high energy asymptotic behaviour of the amplitudes at short distances, 
where one can use the developed methods of perturbative QCD ( pQCD ) to calculate the 
amplitudes. Consider a typical LRG process - the production of two jets with large transverse 
momenta pti ~ — Pt2 ^ yU, with a LRG between the two jets, /x is a typical mass scale of 
"soft" interactions. We have the reaction 

p(l) + p(2) — > Mi[hadrons + jeti{yi,pti)] (1) 

+ LRG[Ay = \yi - y2\] + M2[hadrons + jet2(z/2,Pt2)] 

where yi and 1/2 are the rapidities of the jets and Ay = \yi — y2\ ^ 1- The production of 
two jets with LRG between them can occur because: 

1. A fluctuation in the rapidity distribution of a typical inelastic event. However, the 

Ay 

probability for such a fluctuation is proportional to e~"^, where L denotes the value of 
the correlation length. We can evaluate L «i ^, where ^ is the number of particles 

per unit in rapidity. A LRG means that Ay ^ L and so the probability is small; 



2. The exchange of a colourless state in QCD. This exchange is given by the amplitude 
of the high energy interaction at short distances. We denote it as a "hard" Pomeron 
in Fig.l. We denote by Fg the ratio of the cross section due to this Pomeron exchange, 
to the typical inelastic event cross section generated by gluon exchange ( see Fig. 1 
). In QCD we do not expect this ratio to decrease as a function of the rapidity gap 
Ay = yi —y2- For a BFKL Pomeron [Q, we expect an increase once Ay ^ 1. Using a 
simple QCD model for the Pomeron , namely, that it can be approximated by two gluon 
exchange^, Bjorken[Q] gave the ffist estimate for Fg ^ 0.15, which is unexpectantly 
large. 

As noted by Bjorken, we are not able to measure Fg directly in a LRG experiment. The 
experimentally measured ratio of the number of events with a LRG, to the number of events 
without a LRG ( see Fig. 1 ) is not equal to Fg, but has to be modified by an extra factor 
which we call the survival probability of LRG. 

fgap=<\S\^> xFg. (2) 
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Figure 1: Pictorial definition of fgap, where P and G represent respectively, the exchange of 
a colour singlet and a colour octet. 



The appearance of <| S" p> in Eq. (j^) has a very simple physical interpretation. It is the 
probability that the rapidity gap due to Pomeron exchange, will not be filled by the produced 
particles ( partons and/or hadrons ) from the rescattering of spectator partons ( see Fig. 2a 
), or from the emission of bremsstrahlung gluons from partons taking part in the "hard" 
interaction, or from the "hard" Pomeron ( see Fig. 2b ). 



<| O I > — <| Obremsstrahlung[ ^y — | 2/1 

where s denotes the total cm. energy squared. 



2/2 



> X <| b spectator \ S ) \ > , \6) 



<\ Sbremsstrahiungi Ay ) P> cau bc Calculated in pQCD0, it depends on the kinematics 
of each specific process, and on the value of the LRG . 



To calculate <| Sgpectatoris) p> we need to find the probability that all partons with 
rapidity j/j > yi in the first hadron ( see Fig. 2a ) and all partons with yj < y2 
in the second hadron do not interact inelastically and, hence, do not produce ad- 
ditional hadrons in the LRG. This is a difficult problem, since not only partons at 
short distances contribute to such a calculation, but also partons at long distances for 
which the pQCD approach is not valid. Many attempts have been made to estimate 

's,ec.a.o.(s)P>ii0ii[0] 



<| S. 



11 1, but this problem still awaits a solution. 
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Figure 2: Rescatterings of the spectator partons ( Fig. 2a ) and emission of the 

bremsstrahlung gluons ( Fig. 2b ) that fill the LRG . 



On the other hand, the experimental studies of LRG processes have progressed, and have 
yielded interesting results both at the Tevatron|T2[|T^, and at HERA[]I1|. 



The main results of the experimental data pertaining to the survival probability are: 



the value of <| S p> is rather small. Indeed, the experimental values for fgap{ y/s = 
630 Gel/) = 1.6 ±0.2% ( DOfg) = 2.7 ±0.9% ( CDF|12[ ) and fgapi\^ = 
1800 Gel/) = 0.6 ± 0.2% ( DO[|ll ) = 1.13 ± 0.16% ( CDF[|l| ) can be understood 
only if <| S\^> ^ 1-10%; 



the energy dependence of fgap, namely, the obsevation that 1 12] 1 13 
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leads to <| S" p> which decreases by a factor two in the above range of energy. So, we 
expect that 
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A recent DO estimate 131 of the above ratio is 2.2 ± 0.8. 



It was shown in Ref. ||TT| that the Eikonal Model for the "soft" interaction at high energies, 
is able to describe the features of the experimental data. However, we may question the 
reliability of this approach. Especially worrying, is the energy dependence of <| S p>, since 
a natural parameter related to the parton rescatterings is the ratio 



R 



D 






(4) 



where asD and a^D are the cross sections of single and double diffraction. Experimentally, 
Rd is approximately constant over a wide range of energy. In the Eikonal Model we consider 
o"5D -C a el and aDD -C aei and, therefore, we model Rd — ^ Rei = f^- Experimentally, 
Rei depends on energy, which gives rise to a considerable decrease of the survival probability. 
Indeed, the first attempt to take into account the parameter Rd in the calculation of <| 5* p> 
given in Ref. H] , leads to 



<l^l' 



> 



[1 



2 R 



D 



(5) 



DD = 



Mill 



I I I I I 



Mill 



II II I 



In s 



Figure 3: Lego - plot for double dijfractive dissociation. 



which yields a reasonable value of the survival probability, but cannot account for its sub- 
stantial dependence on energy ( see Ref.|Tl| ). 

The main goal of this paper is to develop a simple model for the "soft" high energy 
interactions, which correctly includes the processes of diffractive dissociation, and to study 
the value and energy dependence of <| Sspectatoris) p>. In the next section we develop 
our approach, while in section 3 we discuss the value of <| Sspectator{s) p> in a model 
which describes the energy dependence of Rei and Rd- We show that our model gives the 
experimentally observed decrease of the survival probability as a function of energy. In 
section 4 we discuss and summarize our results. 



2 A three channel model 



2.1 Diffractive dissociation ( general approach ) 



Diffractive dissociation is the simplest process with a LRG in which no hadrons are produced 
in the central rapidity region. In these processes we have production of one ( single diffraction 
(SD) ) or two groups of hadrons ( double diffraction ( DD ) ) with masses ( Mi and M2 in 
Fig. 3 ) much less than the total energy ( Mi -C ^/s and M2 <^ -^/s ). 



From a theoretical point of view, as was suggested by Feinberg|T^ and Good and 
Walker |1T6[|, diffractive dissociation can be viewed as a typical quantum mechanical process 
which occurs since the hadron states are not diagonal with respect to the strong interaction 
scattering matrix. 



We consider this point in more detail, and denote the wave functions which are diagonal 
with respect to the strong interaction by "^n- The quantum numbers n are called the correct 
degrees of freedom at high energy. The amplitude of the high energy interaction is, therefore, 
given by 

^nma = <'^ni^n2\T\^n[^n'2 > = ^ni.na "^niXi ^n2,n^ ' (^) 

where the brackets denote all needed integrations and T is the scattering matrix. 
The wave function of a hadron is 

oo 

"^hadron = ^^ Cn'^n ■ (7) 

n=l 

For a hadron - hadron interaction we have a wave function "^hadron x "^hadron before the 
collision, while after the collision the scattering matrix T gives a new wave function 

^ final 
J:Z=1 I:;S=1 T.^,=1 T.^,=1 Cn. Cn, < ^n. Vl/„jT|vl/„, Vl/„, > Vl/„, V^,, = 

2^ni=l l^n2=l ^"1 ^1^2 ^ni,n2 ^rii ^ n2 ■ \°) 

Eq. (ID leads to the elastic amplitude 



oo oo 



'^eZ ^ ^ final}^ hadron ^ ^hadron ^ / , / , n\ n2 ni,n2 ; wJ 

ni=l n2=l 

and to another process, namely, to the production of other hadron states, since "^ final may 
be different from "^hadron x "^hadron- 

<yD{s,b) = < "^ finall"^ final > — < "^ finall"^ hadron X "^ hadron > = 

Eoo Y^oo /^2 (^2 /\2 / Y^oo v^oo /^2 (~<2 A \ iMW 

ni=l Z^ra2=l "-"ni "-"m ^711,712 \l^n-i=l Z^ra2=l ^ni ^n2 ^ni,n2 J ■ l-^UJ 

Using the normalization condition for the hadron wave function ( J2n ^n — '^ )•> ^^e cross 
section of the diffractive dissociation processes, Eq. (p!0[) , can be reduced to the form[|l6||jl7| 

aD{s,h) =<\a\s,h)\> " <|a(s,6)|>^ (11) 

where < |/| > = Y.m I]n2 ^ni^n2 fni,n2 ^ud wc havc returned to our original variables: 
energy ( s ) and impact parameter { b ). 



2.2 DifFractive dissociation and shadowing corrections 

For y4„j^„2(s, 6), and only for A„^„2(s,6), we have the unitarity constraint 

2Jm<,„^(.,6) = K,„2(.,6)P + GZ,nMb)- (12) 

Assuming that the amphtude at high energy is predominantly imaginary, we obtain the 
solution of Eq. (|12D 



<,n.(^,fe) = 
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fi - 


^ni,n2(^-b) 
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GZ,n,i^^b) 


= 


1 - 


_ g- f!ni,n2 (»:'') 



; (13) 

(14) 

To find a relation between the processes of diffraction dissociation and the value of the 
shadowing corrections ( SC ), we assume that flni,n2{^^b) ^ 1 and expand Eq. (0) - 
Eq. (p!^ with respect to fini,n2 

K,n2is,b) = ^^^^^^ - %^ + Oinl^^J; (15) 

GZ^,^{s,b) = Qn,,n2{s,b) - ^%M + 0{QI^J. (16) 

Using Eq. (|^) - Eq. (|10D we obtain for the observables 

^ei{s,b) = \ {j:Z=i^Z=iClClQr.,,^,{s,b)y (17) 

O'tot = E^i=l E^=l C'ni Qa ^"i."2('55^) ~ 4 ( 5I]^i=l E^=l C'ni Q2 ^ni,n2 (■^' ^) j (1^) 

Cin = Z^ni=l ^n2=l ^ni ^n2^^™i.'^2 1'^' ") ~ 2 \^ni=l Z^n2=l ^m ^n2 ni,n2 ('^' ") j (1^) 

0"di// = 4(^^=1 E^=l ^^m C'n2^ni,n2('^'^) ^ (5I]^i=l E^=l ^^ni C'n2^"i."2 (■^j &) j }(20) 

In the parton model, one Pomeron exchange corresponds to a typical inelastic event with 
a production of a large number of particles. In this case, we can associate this exchange 
with Q, consequently af^^{s,h) = a^{s,b) oc Q. All terms which are proportional to Q"^ 
describe two Pomeron exchange, and they induce SC. 

We can evaluate the scale of the SC using experimental data on atot, o'ei and <7sd- Indeed, 
we can write the expression for the total cross section in the form 

a,,, = al, - Aaf^' , (21) 

where af^^ is the contribution of Pomeron exchange to the total cross section. Summing 
Eq. (^ and Eq. (^ we derive that 

^^m = ^el + CTdiff = (Tel + (TgD + (^dd = an ■ (22) 
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Figure 4: The experimental data on the ratio R^i = CF^i/atot ( Fig. 4a ) and on the ratio 
Rd = {(^ei + (^SD + (yDD)l(ytot ( Fig. 4b ) versus log(s/so) with Sq = 1 GeV"^. 

The ratio Rd ( see Eq. (^ ) gives the scale of the SC in the situation when these are 
sufficiently weak ( fini,n2 ^ 1 )) since Eq. (|22|) can be rewritten in the form 



R 



D 



C^tot 



(23) 



Fig. 4b shows that over a wide range of energy Rd ~ 0.34 , and it appears to be independent 
of energy. The large value of Rd implies that SC should be taken into account, and that SC 
lead to the small value of the survival probability. The almost constant Rd suggests that 
the SC cannot induce the observed strong energy dependence of the survival probability. 
However, the value of Rd is so large that we have to develop an approach to calculate SC 
for fi„i,„2 ^ 1. 



2.3 The Eikonal Model 



We first discuss the Eikonal model. This is an approximation which has been widely used to 
estimate the value of the SC, in a situation when they are not small. The main assumption of 
this model is that hadrons are the correct degrees of freedom at high energy. In other words. 



we assume that the interaction matrix is diagonal with respect to hadrons. From Eq. (^ 
- Eq. ([np one can see that this model does not include diffractive dissociation processes. 
Therefore, the accuracy of our estimates in the Eikonal Model will be given by the ratio 
(o"sz) + (^DD)/crei- From Fig. 4 one can see that this ratio is about 1 at y/s ~ 20GeV and 
decreases reaching the value of about 0.4 - 0.5 at the Tevatron energies. Thus, we cannot 
expect the Eikonal approach to yield a reasonable estimate. However, this model has the 
advantage of being simple, and it provides a good illustration, given below, of the main 
elements and approximations used in previous calculations. 



1. Assumptions: 

• Hadrons are the correct degrees of freedom at high energies ; 

• (^SD < CTei and aoo < (^ei ] 

• At high energy the scattering amplitude is almost pure imaginary ReUei <^ 
Im aei ; 

• Only the fastest partons can interact with each other . 

The last assumption is the most restrictive, and clearly indicates how far from reality 
the Eikonal Model estimates could be. 

2. Unitarity: 

In the Eikonal Model we only have one amplitude, since the scattering matrix is diag- 
onal in the hadronic wave functions. Therefore, the unitarity constraints of Eq. ([12D 
simplify to 

2Ima,i{s,b) = Ki{s,b)\^ + G'^{s,h) ■ (24) 

Eq. (1^) has the solution 

aei{s,h) = i(^l - e-^) ; (25) 

^^"(5,6) = 1 - e-^(^''') . (26) 

3. The Pomeron hypothesis: 



The main assumption of the Eikonal Model is the identification of the opacity Q{s, b) 
with a single Pomeron exchange, namely 

Q{s,b) = QP{s,b) ; (27) 

nP{s,b) = aoTib) = -^(^)APe-^ = z/(s) e"^ ; (28) 

R\s) = 4Rl + 4a'p\n{s/so) ; (29) 

K^) = ^^ifX" = fi^(^,& = 0). (30) 

We assume a Pomeron trajectory ap(t) = ap(0) + a'pt. Eq. (|27| ) is a reasonable 
approximation in the kinematic region where Q is small, i.e. either at low energies or 
at high energies when b is large. Therefore, the Eikonal approach is the natural gener- 
alization of the single Pomeron exchange satisfying s-channel unitarity. Eq. (pTj) is an 
explicit analytic expression for the well known partonic picture for the Pomeron struc- 
ture, namely, the single Pomeron exchange is responsible for the inelastic production 
of particles which are uniformly distributed in rapidity. 

4. Exponential parameterization: 

In Eq. ( |28| ) a Gaussian form is explicitly assumed for the profile function T{b). This 
corresponds to an exponential form in t space. 

1 



This form is assumed due to its simplicity, the resulting integrals can be done analyt- 
ically and we can write the explicit answer for the physical observables. The proper 
definition of the profile function r(6) is given by the following Fourier transform 

r(6) = ^/ d'b e-^^-"^ GUq') , (32) 

where Gp{q'^) is the Pomeron form factor. For example, one can take for Gp{q'^) the 
prediction of the Additive Quark Model ( see Ref.||18[ ) in which Gp{q^) = Gem, 
where Gem is the electromagnetic form factor of a hadron. 

5. Rd : 

Using Eq. ( pTf ) - Eq. ( pOD one can obtain a closed expression for the "soft" observables 
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atot = 2nR\s){\n{v{s)/2) + C - Ei{-v/2)} ; (33) 

a,n = 'KR\s){\n{v{s)) + C - Ei{-v)} ; (34) 

aei{s) = atotis) - ain{s) ; (35) 

where Ei{x) = /f^ ^ dt and C = 0.5773. 

We define 

P _ ^ ^ Hi\ + C + Ei{-v) - 2Ei{~^) 

'' atot 2[ln[|] + C - E^(-|)] ■ ^""^^ 

Note, that the ratio Rei depends only on z/ and does not depend on the value of radius. 
Using Eq. (|36| ) one can find the value of z/ from the experimental data on Rf.i ( see 
Fig. 4a ) this was done in Ref. |[ll| . With the value of v determined in this way we can 
calculate the survival probability. 

6. Survival probability : 

From Eq. (|26D one can conclude that the factor 

P(s,6) = e-^(^'') , (37) 

is the probability that the two initial hadrons do not interact inelastically. In a QCD 
approach, this means that the fastest parton from one hadron does not interact with the 
fastest parton from another. Therefore, in the Eikonal Model the survival probability 
can be easily calculated in the following way00 

where AHp{^y-, h) is the cross section for a two parton jet production with a LRG 
due to single "hard" Pomeron exchange ( see Fig. 1 ). It has been proven that for a 
"hard" cross section, the h dependence can be factorized out [^ ||19| pO[ • If "^^ assume 
AHp{^yi b) to be Gaussian we have 

AHp{Ay,b) = aHp{Ay)rHib) = ^^^^^ e'^ , (39) 



nR 



where Rjj is the radius of the "hard" interaction. 

Based on this assumption we obtain for the survival probability: 



2^ a7[a,z/] 



<| S spectator{ S ) \ >— , (40) 

where the incomplete gamma function 7(0,0;) = Jq z°'~^e~^dz and a = —^ . 

11 



In Ref . [0 , Eq. (pq) and Eq. (^OD were used to calculate the value and energy dependence 
of the survival probability. It was shown that both the value and the energy dependence 



are sensitive to the value of the "hard" radius, which was extracted in Ref. |11| from the 



experimental data on (i) vector meson diffractive dissociation in DIS at HERAEH], and on 



(ii) the CDF double parton cross section at the Tevatron[^ 



Even though the Eikonal Model, as used in Ref.[|ll[], can reproduce both the experi- 
mentally measured value and its energy behaviour, the reliability of such an approach is 
questionable. In the following we attempt to construct a more realistic model. 



2.4 A three channel model: assumptions and general formulae 

We want to construct a model which takes into account the processes of diffractive dissoci- 
ation, for the case when these are not small ( asn ~ (^ei)- 

2.4.1 Assumptions: 

The main idea behind the model which is presented here, is to replace the many final states 
of the diffractively produced hadrons by one state ( effective hadron ). Doing so, we assume 
that we have two wave functions which are diagonal with respect to the strong interactions: 
\E'i and \E'2. In this case the general Eq. (|^ can be reduced to the form 

^hadron = « ^1 + /? ^2 , (41) 

with the condition a"^ + (3"^ = 1, which follows from the normalization of the wave function. 
The wave function of the diffractively produced bunch of hadrons should be orthogonal to 
"^ hadron and has the form 

■qjjj = -/3^i + a^2 • (42) 

Eq. ( |4^ ) is the explicit form of our assumption, that we replace the complicated final state 
of diffractively produced systems of hadrons by one wave function ^P^j. 



2.4.2 General formulae: 



Substituting Eq. ( |iT| ) and Eq. ( ^ in Eq. (^ we can obtain 

"^ final = a^A.l^l X ^1 + a/3v4i,2[^l X *2 + ^2 X ^i] + /3M2,2^2 X ^2 • (43) 
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The amplitudes Ai^k {i, ^ = 1? 2) can be written in the form of Eq. (|T^) and Eq. (0) ( see 
Fig.5 ) 

At,is,b) = ^{l - e-^M^}; (44) 

Gf^f^{s,b) = 1 - e-^^>fc(^'^) . (45) 

The elastic amplitude is equal to 

(^el\^y 0) ^ ^ hadron ^ ^ hadron I-' | ^ hadron ^ ^ hadron -^ ^ ^ hadron ^ ^ hadron\^ final ^ 

«Mi,i + 2«2^Mi,2 + (3^ A2.2 ■ (46) 

For single diffraction we have the amplitude 

Q'Sois^h) = < "^hadron X "^ hadron \T\ "^ D X ^ hadron > = < "^ hadron X "^ dI"^ final > = 

aP{ -a^Ai^i + («2 _ p^)Ai,2 + /?' ^2,2 } , (47) 

while the amplitude for double diffractive production is 

aDD{s,b) = < "^hadron X "^ hadronlTl"^ D X ^^ > = < ^^ X ^cl^/mai > = (48) 

a^/3^A,^, - 2A,2 + ^2,2} . 



Eq. (^Up - Eq. ( ^5]) together with Eq. ( ^ij) and Eq. (|3) give the general formulae for our 
model. 

In general this, three channel model is an attempt to sum all rescatterings shown in 
Fig. 6. All of them are eikonal type rescatterings, but contrary to the Eikonal Model, the 
quasi-elastic rescatterings with production and rescattering of an effective diffractive state 
has been taken into account. 

We call this model a three channel one because, three physical processes: elastic scat- 
tering, single and double diffraction, are included in it. All general formulae of Eq. (^ - 
Eq. (|TTD were known long ago ( see Ref . 1 2^ , for example ) . In Ref . [§| this general formalism 



was applied to obtain estimates of the value of the survival probability. We now develop a 
systematic analysis to obtain the value of <| S*^ |> in the framework of the three channel 
model, utilizing the experimental data pertaining to the "soft" processes that have been 
measured. 
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1,1/ y^ 

1 1 2 ^2 '2 ^2 

^ ~ 1 1 - exp(-^j j/2) 1 - Qxp(-a^J2) 1 - exp(-a^J2) 



a<< 1 





^,/2 




Figure 5: Solution of the unitarity constraints and the Regge parametrization for the correct 
degrees of freedom in the three channel model. The waved lines denote the exchange of the 
Pomeron. 





Figure 6: The Pomeron interaction in the three channel model. 
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2.5 A three channel model: physical observables 

For the opacities fij ^ we use the same approach as in the Eikonal Model ( see Eq. (|2 
Eq. dD and Fig. 5 ) 

R^ As) _ QkQi ( ^ \ fl^, (s) 



where we have used Reggeon factorization as it is shown in Fig. 5. In this paper, we do not 
use the exact Reggeon dependence on energy, but we utihze the factorization properties which 
appear in Eq. (^). In addition we need to take into account the factorization properties of 
the radii 

Rl, = 2Rl + 2Rl, + Aa'p Hs/s,) , (50) 

where R}q is a radius which describes the t-dependence of the Pomeron - hadron i vertex. It 
is obvious from Eq. (EDf) and Eq. (pDI) that Vti^k can be written in the form 



^1,1 = T^iX- (51) 

1^2.2 = z^2^^ ; (52) 

^1,2 = ^ViV2r{2-r)X'- (53) 



_bf__ 
where X = e ^i.i'''-' and r = ^rVr- 



Eq. ( ^iD - Eq. (|53|) together with Eq. ( ^Sf ) - Eq. (|48|) allow us to express all physical 
observables through the variables i^i, z/2,r and (3. The first three variables depend on energy 
squared ( s ) ( see Eq. (^9]) and Eq. (0) ) while /3 is a constant in our model. From Eq. (^Qf ) 
we expect V2 to be proportional to z/i, and r to be only weakly ( logarithmically ) dependent 
on energy. 

From Eq. (|5T|) - Eq. (|53D we have the following expressions for the amplitudes Aj ^ in 
terms of our variables 

Ai,i(i.i,X) = z (l - e-^) ; (54) 

A2,2{y2,r,X) = ill - e-"-^^^] ■ (55) 



Ai,2(^^i,i^2,r,X) = z 1 - e-^^^. . (56) 
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After some simple calculations we have 

a^^t = 27cRl, J^ ^ {(i-/52)Mi,i + 2(l-/52)/5Mi,2 + P'A2,2} ; (57) 

a^i = nRl,!,' ^{{l-p^A^^, + 2 {1 - (3') ^^ A,^^ + P' A.^^V ; (58) 

asD = 7rRl,{l-(3')(3' J,' ^ {-{1 - p^) A,^^ + {I - 2 f3') A,^, + /5M2,2 F; (59) 

CTDD = nRl,{l-(3^)'(3^ J^ f {Ai,i - 2Ai,2 + ^2,2}' • (60) 

One can see that the ratio Rd as well as Rei, does not depend on Rf^. We will use these 
ratios to fix our variables from the experimental data. 

2.6 A three channel model : survival probability 

To calculate the survival probability of the LRG in our three channel model we recall that 
the physical meaning of the factor 

p^^^ = e-^'.'^^^'^) = {1 - A^,Y (61) 

is, the probability that two hadronic states with quantum numbers i and k scatter, without 
any inelastic interaction at given energy s and impact parameter b. Therefore, we have to 
multiply the "hard" cross section of their interaction by Pi^, and sum over all possible i and 
k for hadron - hadron collisions to obtain the survival probability. The cross section for two 
jet production with large transverse momenta and LRG, can be reduced to the form 



(rH{s,b) = (62) 

[l-PY^i'-^e -M + 2(l-/3^)/3V^i^2a(2-a)5|^e"t + P^u.^e'^-. , 



l^e'^ + 2(l-/3^)/3V-i-2a(2-a)§^e"* + (3'u,§ 
^lA -n^i,2 -n^2,2 



where i?f ^ denotes the "hard" interaction radius of two states i and k. Strictly speaking, 
Rff^ = R'ik{s = So) but, really, we do not know the value of sq- However, we will show in 
the next section that we are able to find the value of Rf,^ directly from experimental data. 

Finally, the survival probability of the LRG is 

<l ^ spectatori ^ ) l> — TTTTV ' ("OJ 



where 



D s) 



rl fix 

N{s) = l^ (64) 
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and 

D{s) = I j^ (65) 

{{l-P^fuiai^iX'^''' + 2(l-/3')z/i,2ai,2X'^^-^'- + /3V2a2,2X'^^-^'/^'""^ } , 



where we denote 1^12 = J ui 1^2 r {2 ~ r) and ai^k = -p^- In the next section we will de- 

termine all parameters from the experimental data, and will find the value and the energy 
dependence which are typical for the survival probability in our model. 

3 Numerical evaluation of <| Sg ej^^^tor l> from the ex- 
perimental data 

3.1 Fixing the parameters of the model 



Following the ideas of Ref. |Tl| we determine all the parameters of the model directly from 
the experimental data. 

1. The most striking experimental fact is that Rd is almost independent on energy. Rd 
is rather big ( Rn ~ 0.4 ) ( see Fig. 4b ). We found that this energy behaviour, as 
well as the value of Rd, allowed us to find the coefficient t in the equation P2 = ti^i. 

2. The energy behaviour and the value of Rei ( see Fig. 4a ) is used to determine the 
value of z/2 at different energies, as has been done in Ref. [ 



3. Unfortunately, no reliable measurement has been performed on the double diffractive 
cross section. We used the estimates of Ref. [^ for the value of this cross section, to 
check that our choice of ui and z/2 is not in a contradiction with the value of the double 
diffraction cross section. 

4. We used the experimental data on hard diffraction in DIS to fix the value of the "hard" 
radii in Eq. (|6^). 

5. Our goal is not to fix all parameters, but rather to find out how sensitive the value and 
the energy behaviour of the survival probability is, to uncertainties in the values of the 
model parameters. We also evaluate the range of typical values of <| S^p^^^^^^^^. >. 
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Fie;. 7-a 



FiK.7-b 



Figure 7: Ratios Rei = o'ei/o'tot and Rd = {cei + <Jsd + o'DD)/o'tot versus vi and t, where t 
is defined as 1^2 = t ^i- 

3.2 Rd and Rei 

We start by fixing the model parameters for the case of very high energies for which r — > 1. 
In Fig. 7 we plot Rd and Rei at fixed /3 = 0.65 versus ui and t, where t is defined z/2 = tui. 
We have argued that in the Regge approach we expect that 1^2 is proportional to z/i. One 
can see from Fig. 7, that Rd is about 0.4, only for large values of t. Note that z/i 2 = Vii^i 
at high energies. Therefore, for t > 100 we have 1/12 > lOz/i in accordance with the global 



fit of the experimental data on "soft" processes [P4 



In Fig. 8 we take t = 300 and plot Rd and Rei versus ui . One can see that Rd is a 
very smooth function,, while Rei depends substantially on z/i. Such a behaviour refiects the 
experimental situation shown in Fig. 4. We use the experimental data for Rei given in Fig. 4 
to assign a definite value of z/i for a definite value of energy. In particular, we find ui = 0.25 
for v^ = QAOGeV and z/i = 0.5 for ^ = 1800 Gel/ to give Rei = 0.187 and Rei = 0.237, 
respectively. 

Comparing Fig.8a and Fig.8b we can see that /5 - dependence for P > 0.5 is not essential. 



3.3 "Hard" radii 



Before calculating the value of the survival probability we discuss the values of the "hard" 
radii R^j^. Fortunately, these can be determined directly from the experimental data on 



diffractive production of vector mesons in DIS|2^ at HERA. The data show that this pro- 
duction depends differently on the momentum transfer for elastic ( see, for example, the 




Fie. 8-a 



Fig.8-b 



Figure 8: Ratios R^i = cTei/cftot o,nd Rd = (cei + ctsd + o'DD)/o'tot versus z/i for two values 
of 13 ( [3 = 0.65 ( Fzg.Sa) and p = 0.5 ( Fig.Sh ). 



process 7* + p — > ^ + p in Fig. 9 ) and inelastic ( 7* + p -^ \I' + X in Fig. 9 ) production. 
In the exponential parameterization, the t - dependence are characterized by the slopes Bf.i 
and Bin. 

Experimentally, these slopes are B^i = 4Gel^~^ and Bin = l.QQGeV~^. Since the 
vertex 7* ^ \E' does not depend on t at large value of photon virtuality Q^, we can view 
such an experiment as a way of measuring the t - dependence of the Pomeron - hadron form 
factor, and/or the transition form factor of a hadron to a diffractive state, due to Pomeron 
exchange. Incorporating the experimental data on slopes in the expressions of the "hard" 
radii we find 



^1,1 



R? 



2 

I 
'2,2 



16 GeV-^ ; 
11.32 GeV-^ 
6.64 GeV-^ . 



(66) 



One can see from Eq. (^) the value of the survival probability depends on the ratios rj fc. 
To calculate these ratios we have to specify the values of "soft" radii R^^.. We assumed the 
following values for the "soft" radii 



-^1,1 



12 + 4a^(0) ln(s/so) GeV- 
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Figure 9: Two slopes in diffractive J/\I/ production in DIS in the Additive Quark Model. 



Rl^ = 6 + 4a^(0) ln(s/so) Ge\/-2 
Rl^ = 4a'p\n{s/so) GeV-K 



(67) 



It should be stressed that these radii are taken from our parameterization of the "soft" data, 
but they also describe the data on the elastic slope, and agree with known information on 
the t - dependence of single diffraction dissociation||23l . 



3.4 <| Ogpgptatorl^) l> 

In Fig. 10 ( solid curve ) we show the value of the survival probability at ^/s = 1800 GeV as 
a function of (3. This value turns out to be rather small, in agreement with the experimental 
data Ik 



Fig. 11 ( solid curve ) shows the energy dependence of <| Slp^^^-i-^[ 
ratio 

spectator \ 



R 



<\ SL^,.,„,„A ^fs = 640 GeV) \^> 



^1 c2 ( r: 

^ I ^ spectator v V "^ 



1800 Gel/) P> 



s) p>, namely the 



(68) 



versus p. 
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In both these figures, at each fixed /3, we found the values of z/i for i/i = 640 GeV and 
for ^/s = 1800 GeV^ from the values of Rei ( Rei = 0.187 and Rei = 0.237, respectively, see 
Fig. 4a ). The value of the survival probability was calculated using Eq. (|6^). 

One can see that the ratio Rs reaches the value oi 2.2 for (3 = 0.6, which we consider a 
typical value which fits the experimental data. It should be stressed that for /? < 0.5 the 
value of Rd is smaller than 0.3. This fact is in contradiction with the experimental data 
shown in Fig.4b. At /5 = we reduce to the usual Eikonal Model and we obtain a ratio Rs 



which is about 1.6. In Ref.[ll] we found a larger value as we used an average value for the 



"hard" radius, while here we have introduced different radii for the different processes. 

We would like to stress, that the accuracy of our estimates is not high. This is mostly 
because of the large dispersion of the experimental data for Rei ( see Fig. 4a ). To illustrate 
this point we plot the value of <| S*^ |> in Fig. 10 ( dotted line ) taking Rei = 0.215 at 
a/s = 1800 GeV. One can see that the difference between these two curve is about a factor 
of 2. The spread of the experimental data influences dramatically the ratio of Eq. (|68[) ( see 
Fig. 11 ). To illustrate it we plot in Fig. 11 several lines that correspond to different choices 
oiR. 



el 



Rei{y^ = 640 Gel/) = 0.187 and Rei{^ = 1800^61/) = 0.237 ; (69) 

Rei{y^ = 64:0 GeV) = 0.212 and Rei{^s = 1800 GeV) = 0.237 ; (70) 

Rei{^ = 640 GeV) = 0.187 and Rei{^s = 1800 GeV) = 0.215 . (71) 



The full line in Fig. 11 corresponds to parameters of Eq. (|69D, the dashed line to parameters 
of Eq. (^) and the dash-dotted line to parameters of Eq. ([71|). 

We conclude, therefore, that we cannot obtain a definite prediction even within the 
framework of a particular model, due to the uncertainties in the experimental data on Rei. 



4 Discussion and Summary 

In this paper we give an example of a model in which 



1. The processes of elastic and diffractive rescatterings have been taken into account, 
their cross sections being of the same order of magnitude ( aso ~ cTei and ctdd ^ aei 

); 
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Figure 10: The value of <\ S"^ \> at y/s = W = ISOOGeV^ versus (3 in the three chan- 
nel model. The full and dashed lines corresponds to parameters of Eq. ( [g^ j and Eq. (fT^j, 
respectively. 
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Figure 11: The value of ratio <\ S'^ |> {W = 640 GeV)/ <\ S^ \> {W = 1800 Ge\/) versus 
(3 in the three channel model. W = ^/s. The full, dashed and dash-dotted lines corresponds 



to parameters of Eq. (WWj Eq. (iTOj and Eq. (71), respectively 
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2. It was shown that the scale of the SC is not given by the ratio Rd, but rather by 
the separate ratios aei/cJtoti <^SDl<^tot and OBDJatot- Since each of these ratios shows 
considerable energy dependence, we do not expect a constant survival probability, 
contrary to the simpler model of Ref. |0; 

3. It was demonstrated that the small value of the survival probability, as well as its 
strong energy dependence appear naturally in our approach; 

4. The rather large value of V2 ~ 300 v\ reflects (i) a smaller value of W^ 2 i^ comparison 
with i?f ^ observed experimentally, and (ii) the fact that this value takes into account 
the integration over the mass of the produced hadrons in our oversimplified model; 

5. The parameters that have been used are in agreement with the more detailed fit of the 
experimental data on "soft" processes ( see Ref.[|^ ). 

Theoretical predictions for the value of the survival probability are still not very reliable. 
However, developing different models enables us to learn and assess which class of models 
provides natural predictions for both the value and the energy behaviour of the survival 
probability. We want to draw the readers attention to the fact that our estimates for the 
value of the survival probability given in Figs. 10 and 11 are very close to the estimates 
obtained in the Eikonal Model |jll|, in spite of the fact that the three channel model is quite 
different from the eikonal one. 

New measurements both on LRG processes, and on the cross sections of diffraction disso- 
ciation ( in particular, on double diffraction ) would be very useful for a deeper understanding 
of "soft" interactions at high energy. 
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